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Decay of the monochromatic capillary wave. 
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It was demonstrated by direct numerical simulation that, in the case of weakly nonlinear capillary waves, 
one can get resonant waves interaction on the discrete grid when resonant conditions are never fulfilled exactly. 
The waves's decay pattern was obtained. The influence of the mismatch of resonant condition was studied as 
well. 
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Nonlinear waves on the surface of a fluid are one of 
the most well known and complex phenomena in na- 
ture. Mature ocean waves and ripples on the surface 
of the tea in a pot, for example, can be described by 
■ ^ "very similar equations. Both these phenomena are sub- 
^~^iStantially nonlinear, but the wave amplitude is usually 
significantly less than the wavelength. Under this con- 
I ^^ 'dition. waves are weakly nonlinear. 

I To describe processes of this kind, the weak tur- 
'bulence theory was proposed It results in Kol- 

.mogorov spectra as an exact solution of the Hasselman- 
Zakharov kinetic equation [31. Many experimental re- 
sults are in great accordance with this theory. In the 
case of gravity surface waves, the first confirmation was 
■obtained by Toba and the most recent data by 
|Hwang were obtained as a result of lidar scanning 
■of the ocean surface. Recent experiments with capillary 
iwaves on the surface of liquid hydrogen 0,171 are also 
'in good agreement with this theory. On the other hand, 
.some numerical calculations have been made to check 
'the validity of the weak turbulent theory [5).p|.pi1]. 
In this Letter we study the one of the keystones of 
. 5^ the weak turbulent theory, the resonant interaction of 
weakly nonlinear waves. The question under study is 

H the following: 
^ , 

• How does a discrete grid for wavenumbers in nu- 
merical simulations affects the resonant interac- 
tion? 
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We study this problem for nonlinear capillary waves 
on the surface of an infinite depth incompressible ideal 
fluid. Direct numerical simulation can make the situa- 
tion clear. 

Let us consider the irrotational flow of an ideal in- 
compressible fluid of infinite depth. For the sake of sim- 
plicity, let us suppose fluid density p — 1- The velocity 
potential (f) satisfies the Laplace equation 



A0 = 

in the fluid region bounded by 
—oo < z < ri{r), : 
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with the boundary conditions for the velocity potential 
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on z = ?7, and 
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on z ^ — oo. Here r] = rj{x, y, t) is the surface displace- 
ment. In the case of capillary waves, the Hamiltonian 
has the form 

H = T + U, 



• Can a nonlinear frequency shift broad resonant 
manifold to make discreteness unimportant? 
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(fr / {V(t>fdz, 



(5) 
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U = (7 (v/l + (Vr;)2 - l)(fr, 



(6) 



where cr - is the surface tension coefficient. In 

it was shown that this system is Hamiltonian. The 
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Hamiltonian variables are the displacement of the sur- 
face ri[x,y,t) and velocity potential on the surface of 
the fluid tp{x,y;t) = (j){x,y,ri{x,y]t)]t). Hamiltonian 
equations are 

d-q _5H d'4) _ SH 
dt Sip ' dt Si] 

Using the weak nonlinearity assumption 3 one can ex- 
pand the Hamiltonian in the powers of surface displace- 
ment 
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The third order is enough for three-wave interactions. 
Here, k is the linear operator corresponding to multi- 
plication of Fourier harmonics by the modulus of the 
wavenumber k. Using Q, one can get the following 
system of dynamical equations: 



f] — kip — div(?7VV') — k[f]kil>\, 
1 r 



ij) — aArj - 



(9) 



The properties of fc-operator suggest exploiting the 
equations in Fourier space for Fourier components of 
rj and ip, 
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Let us introduce the canonical variables at as shown 
below 
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where 

LJk - V^. (11) 

With these variables, the Hamiltonian (jH)) acquires the 
form 



H = ujk\ak\^dk+ 
11/ 

k2 + ko)dkirfk2rfko + 



X(5(ki 
1 1 
"^22^ 
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^^k°k2 ("^ki ak2 flko + «ki «ko ) X 

k2 - ko)dkidk2dko. 



Here, 
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-^kik2 = (kik2) -I- |A:i||A;2| 



The dynamic equations in this variables can be easily 
obtained by variation of Hamiltonian 

5H 

6al 
i 1 f 

/ ^kik2"kiak2'5(ki + k2 - k)dkidk2- 

z / k ^^^^ 
/ M^^^al^akJ{k + k2 - ko)dk2(iko- 

/ ^^ik2«2^(ki+k2+k)dkidk2. 

Each term in this equation has its own clear physical 
meaning. The linear term gives a periodic evolution of 
the initial wave. The first nonlinear term describes a 
merging of two waves ki and k2 in k. The second de- 
scribes a decay of the wave kg to the waves k and k2. 
And the last term corresponds to the second harmonic 
generation process. It is useful to eliminate the linear 
term with the substitution 

Ok - Ai^e'^"*. (15) 

In this variables, the dynamical equations take the form 



2 27r . 

X(5(ki + k2 - k)dkidk2- 
-h J <k2^k2^koe-^<^ 
X(5(k -I- k2 - ko)(ik2dko, 



(16) 



where 



^felfc2 ~ + ^k2 - ^ko- 



(17) 



Here we do not consider the harmonic generation term. 
The remaining terms give us the following conditions of 
resonance 



^fci + ^^fe2 - ^/c = 0, ki-Hk2-k = 0. (18) 



All this theory is well known in the literature 

Now let us turn to the discrete grid. Also, from this 
point we assume periodic boundary conditions in x and 
y with lengths Lx and Ly. One can easily obtain equa- 
tions similar to (|16|l 



kik2 



xA 



(ki+k2) ,-k- 



(19) 



i27r 

L'xLy 



^'■^kk2"^k2^koe 2 ^(k+k2).-ko: 



k2ko 



where Aki,k2 is the Kronecker delta - the discrete ana- 
logue of the Dirac delta function. 
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Consider the decay of a monochromatic capiUary 
wave ylk^ on two waves 
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Let Aki,Ak2 be small (|ylko| > max(|^ki |, l^ka I) at 
t = 0). In this case the equations can be linearized. The 
solution of linearized l|2()|l has the form (Ak,, ~ const) 



where 



^ki,.(i) = Aki.,(0)e 



2 fcife 



At 



(21) 



27r 



I 2 fei*:2 



(22) 



In the case of a continuous media, resonant condi- 
tions (|18|l can be satisfied exactly. But on the grid, there 
is always a frequency mismatch rife^j.^ 7^ although if 
the amplitude of the initial wave is high enough there 
are resonances even on a discrete grid. But the width 
of this resonance is very important. 

System of equations can be solved numerically. 
This system is nonlocal in coordinate space due to the 
presence of the fc-operator. The origin of this opera- 
tor gives us a hint to solve JHJ in wavenumbers space 
(i^-space). In this case we can effectively use the fast 
Fourier transform algorithm. Omitting the details of 
this numerical scheme, we reproduce only the final re- 
sults of calculations. 

We have solved system of equations (O numerically 
in the periodic domain 27r x 27r (the wave-numbers kx 
and ky are integer numbers in this case). The size of 
the grid was chosen as 512 x 512 points. We have also 
included damping for waves with large wave numbers. 
In X-space damping terms for 77k and T/^k respectively 
were the following: 7kfyk and 7k^k, where 7k was of the 
form 



1, 



7k 



7k = 



0,|k| < -|k 
1 ^ 



kmaxl) , |k| > r)|kmaa;|; 



(23) 



here, 70 is some constant. 

As an initial conditions we used one monochromatic 
wave of sufficiently large amplitude with wave numbers 
ko {k[jx = O,koy — 68). Along with that there was a 
small random noise in all other harmonics. 



Resonant manifold (|18|1 for decaying waves 
ko = 

^ ko - ky j ' 

is given at Fig^ Since the wave numbers are integers 





(24) 
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Fig.l. The resonant manifold for ko — 68. 

the resonant curve never coincides with grid points ex- 
actly. A detailed picture is given in FigI21 It is clear 
that some points are closer to the resonant manifold 
than others. This difference might be important in nu- 
merics. 
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Fig. 2. Different mismatch is seen at different grid points. 

In the beginning, one can observe exponential 
growth of resonant harmonics in accordance with (|21(l 
and 1)22(1 . This is shown in Fig|31and Fig01 Here one 
can clearly see that some harmonics are in resonance 
and others are not. 



4 



A. I. Dyachenko, A. O. Korotkevich, V. E. Zakharov 



1 1 1 1^1 1 

kn = (00, 68) 




k = (26,34) 




it = (26, 29) 




k = (26,28) —- 




k = (26, 20) 


1 1 


1 1 1 1 



3 4 
Time 




Fig. 3. Evolution of various harmonics for decaying wave 
ko = (00,68). 



Than almost all harmonics in the resonant manifold 
become involved in the decay process (FigEJ. Later, the 




Fig. 5. Secondary decays start. Time t=ll. 
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Fig. 6. The level lines for |afcp. Secondary decays are 
clearly seen. Time t=14. 



Fig. 4. Resonant harmonics starting to grow. Time t=1.4. 

harmonics that are the closest to the resonant manifold 
(compare with Fig|2Il reach the maximum level, while 
the secondary decay process develops. Waves ampli- 
tudes became significantly different. The largest ampli- 
tudes are for those waves with the maximal growth rate. 
One can see the regular structure generated by the ko 
wave in Fig|H| After a while the whole fc-space is filled 
by decaying waves, as shown in Fig[3 

Direct numerical simulation has demonstrated that 
the finite width of the resonance makes discrete grid 
very similar to continuous. Of course, this is true only 
if the amplitude of the wave is large enough, so that 
according to 



2tt 



> 



1 



(25) 



As regards numerical simulation of the turbulence, 
namely, weak turbulence, the condition (|25|l is very im- 
portant, has to be treated as the level of turbulence. 
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Fig.7. Wave numbers spectrum at time t=57. 
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